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CMOS-hybrid arrays have become competitive optical detectors for use 
in ground- and space-based astronomy. Inter-pixel capacitance is one source 
of error that appears in most CMOS arrays. In this paper we use a single 
pixel reset method to model inter-pixel capacitance (IPC). We combine 
this IPC model with a model for charge diffusion to estimate the total 
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1. Introduction 

Complementary Metal-Oxide Semiconductor (CMOS) sensors have become a competitive 
astronomical ground- and space-based detector solution. CMOS sensors have a flexible read- 
out structure that allows a single pixel or group of pixels on the array to be read out or reset 
at any time without disturbing or reading out the rest of the array. This random-access, non- 
destructive read capability is ideal for dynamic range-driven astronomical applications, since 
it permits bright and faint objects to be observed simultaneously using a single detector. In 
addition to the flexible readout, CMOS sensors are naturally less sensitive to radiation than 
more traditional detectors like Charge Coupled Devices (CCDs), since damage to one pixel 
in the array does not adversely affect subsequent pixels in a row or column of the array. This 
inherent radiation hardness is particularly appealing for space-based applications. 

While the readout capabilities of CMOS sensors are ideal for a variety of observing strate- 
gies, the flll factor of each pixel is signiflcantly reduced because the readout circuitry is 
implanted directly on the photodetector material. Janesick et al. [1] describe how the CMOS- 
hybrid focal planes have been developed to address this issue. A CMOS-hybrid sensor is a 
CMOS device or Readout Integrated Circuit (ROIC), mated with a layer of photodetector 
material. The two layers are typically joined together using indium bump bonds. The re- 
sultant hybrid SCA is back-illuminated, and combines the flexible readout of the CMOS 
with a CCD-like flll factor of 100%. The United States Naval Observatory (USNO) has 
used CMOS-hybrid detectors to take advantage of this performance and flexibility, as part 
of its development of very large format focal plane technologies, in support of the Joint 
Milli-Arcsecond Pathflnder Survey (JMAPS) astrometry mission. 

USNO has been testing large format, Teledyne Imaging Sensors (TIS) H4RG-10 Hybrid 
Visible Silicon Imager (HyViSI) Sensor Chip Assemblies (SCAs) since the development of 
the flrst generation-Al detector in 2006, as described by Borland et al. in 2007 [2]. In 2008 
USNO supported development of a second generation- A2 H4RG-10 with signiflcantly lower 
dark current than the Al predecessor, again described by Borland et al. in 2009 [3] A third 
generation detector, the H4RG-10 BO detector, was fabricated in 2010 with JMAPS support 
to address yield and pixel operability requirements. This third generation of detector has very 
low dark current and noise properties comparable to many CCBs. The quantum efficiency 
is in excess of 80% across multiple wavelengths and the non-linearity is < 1%. Additionally, 
with JMAPS risk reduction support, Hubbs et al. in 2011 [4], describes how TIS was able 
to increase pixel operability to better than 99.9%. Here, operabihty specifically refers to the 
percent of pixels that are fully connected without shorting. 

While these fiexible and low-noise CMOS hybrid arrays are excellent for most astronomical 
applications, the H4RG-10 show higher levels of crosstalk, or more specifically: the Inter-Pixel 
Capacitance (IPC) component of crosstalk than CCBs. We also note that some CMOS pixel 



circuit designs do not exhibit high levels of IPC. However, the chosen design of the H4RG-10 
has advantages that are difficult to achieve with these other designs, namely low read noise, 
low dark current and low power consumption. Indeed, most low background applications 
prefer the source-follower approach despite the IPC problem (eg. JMAPS H4RG-10 BO), 
while higher background applications typically use one of the other design options. 

As discussed in detail below, crosstalk can be problematic for applications like astrometry, 
because it has the effect of blurring the point spread function (PSF) of the photon source, 
resulting in lower effective signal to noise (S/N) in each pixel and high centroiding errors 
(See also [SHT]). For this reason JMAPS has supported development of a fourth generation 
detector, the H4RG-10 Bl, designed in part to reduce IPC to values that are negligible for 
astrometry. 

In this paper we discuss the data analysis and modeling that have been used to understand 
crosstalk for astronomical detectors. In Section 2 we describe the models for charge diffusion 
and IPC that were combined to create a ^^Fe model based on the H4RG-10 BO detector 
(third generation). In Section 3 we describe the JPL camera used to collect the data. In 
Section 4 we discuss the analysis methodology for single-pixel-reset data and compare the 
results with the modeled data. In Section 5 we discuss the analysis methodology used for ^^Fe 
data obtained with the astrometric camera, and again compare the analysis and modeling 
results. Finally, in Section 6 we summarize our findings. (In a subsequent paper we show 
how crosstalk affects astrometry and photometry using detailed detector and optics models 
for a realistic astrometric telescope (See [5] for details.)) 

2. Crosstalk models 

2. A. Introduction to charge diffusion and IPC 

Crosstalk between pixels is caused by two independent phenomena, charge diffusion, and IPC. 
Charge diffusion is the lateral movement (i.e. pixel-to-pixel) of charge between the points 
of charge production and charge collection in the bulk substrate of the detector. Charge 
diffusion occurs in all photosensitive material, including the silicon substrate of all CCDs. 
IPC is the capacitance that arises between adjacent detector pixels in the source-follower 
CMOS design, and leads to coupling of signal between those pixels via displacement currents 
flowing from the collection node. 

For optical detectors with photosensitive silicon material, crosstalk can be measured by 
exposing the detector to a radioactive source with a known charge production rate in bulk 
silicon, and then measuring the voltage in pixels surrounding a central hit. ^^Fe is a soft X- 
ray source, quickly decaying (half-life: 2.7 years) to Mn when a K shell electron is absorbed 
into the nucleus. An electron from either the L or M shells drops to flll the hole created in 
the K shell. This drop to a lower energy level causes the emission of either a Kq, (5.9 KeV) 



or K^ (6.5 KeV) X-ray. The absorption of a K^ photon in the bulk sihcon of the detector 
produces, on average, 1620 e~, while K^ photon absorption produces an average of 1778 e~. 

In ^^Fe testing, the detector, under controlled environmental conditions, is briefly exposed 
to the radioactive source and the resulting charge is collected and read out for each exposure. 
When analyzed, photon hits approximately normal to the detector surface and centered on a 
single pixel are averaged and normalized to create a kernel representative of pixel crosstalk. 
^^Fe testing and analysis is described in more detail in Section 5 of this document. 

In addition to ^^Fe testing, single-pixel-reset testing (SPR) allows for the direct character- 
ization of IPC alone, requiring no illumination source [7]. In SPR, after setting all pixels in 
the SCA to a single voltage and making an initial readout, a well-spaced grid of single pixels 
is reset to a second voltage level. A subsequent readout of the pixels will reveal any IPC as 
signal in pixels adjacent to the reset pixels. As with ^^Fe testing, multiple SPR test results 
are averaged (discarding any bad pixels, edge-affected pixels, etc.) to create a representative 
detector IPC kernel. SPR testing and data analysis methodology is described in more detail 
in Section 4. 

The IPC kernel obtained through SPR testing, when combined with charge diffusion mod- 
eling, can be used for ^^Fe test verification and prediction. This resultant ^^Fe model is 
based on the convolution of a charge diffusion kernel, resulting from an incident Kq, photon, 
with an IPC kernel built on SPR-measured IPC. This new, modeled kernel can be used to 
approximate the average pixel crosstalk expected in ^^Fe testing for a given detector. The 
subsections below describe the modeling of IPC and charge diffusion. 

2.B. IPC model 

A simple model of the IPC expected within a detector array can be developed based on the 
approach described by Moore, et al. in 2004 [8]. The detector array of photodiodes is modeled 
as an array of capacitors, as shown in Figure [1] , each identical and with a capacitance that 
is unchanging with voltage. 

The detector is constructed in a way that allows an electrical field to exist between neigh- 
boring collection nodes, essentially creating small coupling capacitors between the nodes. 
Charge entering a single nodal capacitor Qtotai causes a voltage change in that node, and 
through the coupling capacitors, causes voltage changes in neighboring nodes. Qtotai is the 
sum of all the apparent charge seen both in the voltage of the original node and the voltages 
of n neighboring nodes. 

T.Vn = ^ (1) 

In Moores approach, the impulse response of each node h{n) is a ratio of the charge that 
appears electrically in a node Qn, to the photocurrent that entered the original node Qtotai, 



hin) = -^. (2) 

Using this approach, and assuming two-dimensional symmetry in nearest neighboring pixels 
and diagonal neighboring pixels, a simple model of IPC can be constructed. Charge appearing 
electrically in surrounding pixels is defined in terms of a single variable a, defined as the 
percent of total charge seen in any of the four nearest neighbor pixels. Symmetry in nearest 
neighbors and diagonal neighbors is assumed in this model. An example is shown in Figure 
[2] for a 3x3 pixel array. 

The sum of the charge in all pixels in the array = 1.0. The above model kernel, with a = 
0.075, is an IPC model for a general CMOS detector. The model can be extended to larger 
kernel models to account for detectors with IPC that is more broadly spread. 

2.C. Charge diffusion 

Lateral charge diffusion occurs while charge moves between the point of generation and the 
point of collection, in the detector substrate. The process begins with the absorption of a 
Kq, photon, which, for ^^Fe, produces a cloud of 1620 charge pairs within the substrate, as 
illustrated in Figure [31 

The photon attenuation length la for a K^ (5.9 KeV) photon is: 

^"^ " 47rlm(n)' ^^^ 

using the absorptive term Im(?T,) of the complex refractive index of silicon, n 

n = l-5-il3, (4) 

where 6 is the refractive index decrement and (3 is the absorptive index. For a photon with 
£■=5.9 KeV, /3 ~ 6.0 x 10~^ [12] and la ~ 28/im. In the ^^Fe model, charge diffusion is 
calculated for thin slices of the bulk Si substrate. The probability of photon absorption is 
calculated for each slice: 

Paz = ]- ( e 'a - e 'a j (5) 

where di and ^2 are positions along the z axis (a direction normal to the pixel array plane 
(xy plane), as defined for illustration in Figure [3] above), and di is further from the detector. 
(In this model, the incident photon is always normal to the detector surface. This orientation 
was selected to mimic the symmetry selection criteria for good hits used in actual ^^Fe test 
data analysis, as described below in Section [4.01 ) 



We now consider charge diffusion for a detector with a fully depleted substrate. Charge 
diffusion, as a function of depth and weighted for absorption probability, is well represented 
by the two dimensional Gaussian 
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with a (z) the Root Mean Square (RMS) standard deviation of charge spreading in x and y. 
This spreading is described in terms of the diffusion constant Dp and the transit time ttr for 
charges (holes, for the detectors under consideration) to move from the point of charge pair 
generation to the point of collection 



a(z] = \/2Dpttr. (7) 

The transit time ttr can be found using the hole drift velocity f drift as described by [9|[T0]: 
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and integrating over the entire depletion depth. 
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where esi is the dielectric constant for silicon (11.9 x 8.854 x 10~^^ C/V), q is the fundamental 
charge unit (C), and A^'d is the doping density for the bulk silicon substrate. -Emax includes 
the effects of both the substrate bias voltage V, and the depletion voltage Vd [TT] . 
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where the Einstein relation, Dp/fi = ksT/q has been applied to equation [71 above, and 
where T is the temperature of the detector (K) and /cb is Boltzmanns constant (J/K). For 
absorption of the photon at a specific z within the substrate: 
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In the model, the charge diffusion function /cd is calculated for each incremental slice through 
z, using the above equation for a (z), with: T = 193 K, V = 40 V, and zj) = lOO/xm. 

The charge cloud produced by absorption of a K^ photon is resolved at a higher-than-pixel 
resolution, and diffusion to each neighboring pixel is determined by the position of the cloud 
center projected somewhere on the surface of the central pixel. Because of this, the charge 
diffusion is calculated at a subpixel resolution, with an absorption probabihty weighting for 
each shce. The high-resolution kernel is then summed and rebinned to a lower resolution 
kernel, ready for convolution with the IPC kernel. Examples of the charge diffusion model 
to both subpixel and pixel resolutions are shown in Figure HJ 

3. JMAPS test camera 

The data described here were taken at the Jet Propulsion Laboratory (JPL) using a test 
camera specifically designed for ground based astrometric testing of H4RG-10 detectors for 
JMAPS. The camera functions with 32 outputs controlled by a TIS non-cryogenic SIDECAR 
(system image, digitizing, enhancing, controlling, and retrieving) ASIC (Application Specific 
Integrated Circuit). The measurements were taken at 193K and the detector substrate voltage 
was set to 10, 20, 30 or 40V depending on the measurement. The camera voltage settings 
were the same for IPC SPR and ^^Fe data acquisition. Table [U lists the primary voltage 
parameters used for these measurements. An image of this camera is shown in Figure El 



Table 1. JMAPS Camera Voltage Parameters 



Parameter 



Symbol 



Vohage (V) 



Bias Voltage 

Digital Positive Power Supply 

Analog Positive Power Supply 

Drain Node of Pixel SF 

Drain Node of Output SF 

Source Node of Internal Current for SF 

Detector Substrate Voltage 

Detector Reset Voltage 

Bias Voltage for Current Source of SF 



VSUB 


10-40 


VDD 


3.12 


VDDA 


3.12 


CELLDRAIN 


0.00 


DRAIN 


0.00 


ABIASPOWER 


3.00 


DSUB 


0.50 


VRESET 


0.32 


VBIASGATE 


2.00 



4. Single pixel reset (SPR) data 

4- A. SPR data acquisition and analysis 

The JMAPS astroinetric test camera is equipped witli microcode that can perform the single 
pixel reset (SPR) technique described by Seshadri [12] and Finger [13]. The single pixel reset 
function enables reset of single pixels in the array to a different reset voltage than the rest 
of the pixels in the array. The effective result is charge on the integration node for these 
reset pixels. This method is useful for isolating IPC, since the charge is not generated in the 
photodetector material, and is therefore not susceptible to the effects of charge diffusion. 

IPC maps are created using the following data acquisition method. First the entire detector 
is reset to a given reset voltage and the pixels are read out to generate an offset frame. Then, 
a grid of widely spaced single pixels is reset to a new voltage. Following this reset, all pixels 
are again read out. If IPC is present, the difference of these two images will show voltage in 
pixels neighboring the reset single pixels. A cropped single pixel reset image of a section of 
an H4RG-10 BO detector taken with the JMAPS camera is shown in Figure El 

For this analysis the initial reset voltage was first set to 300mV for all pixels. The SPR grid 
of pixels was then reset to 0, 100, 200, 300, 400, 500, 600mV. This resulted in 7 measurements 
of varying signal levels per substrate voltage (lOV, 20V, 30V, 40V). To analyze the data, the 
following procedure for finding good hits was adopted: 

1. First the image was offset (bias) corrected. 

2. A small range of central pixel values were identified as hits. The range of values for the 
central pixel was identified as a Gaussian distribution in the histogram of the image, 
and was well above the noise floor of the image. Defining this range is important to 
ensure that hot or bad pixels are excluded from the IPC estimate. On average, each 
image contains approximately 64,000 good hits that are used for the analysis. 

3. All pixels defined as hits based on the criteria above were background corrected using 
a local background correction method. An annulus of 2 pixels outside the central hit 
was extracted. The median of the pixels in this annulus was subtracted from the values 
of the inner 9x9 pixels. An image illustrating the annulus and kernel is shown in Figure 

El 

4. The background-corrected kernels were then stacked and averaged. 

5. The resultant measured IPC kernel was normalized, so that the charge in a single pixel 
was represented as a percent of the total charge in the kernel. 

For this detector, the noise floor was reached in both IPC and ^^Fe crosstalk data at the 9x9 
pixel boundary. 



4.B. IPC data and model 

4.B.I. Data 

Figure [8] shows the fraction of signal in the central pixel as a function of the SPR voltage. 
In this plot, the saturation limit is indicated by the drop-off in fractional signal at 400 mV. 
However, in order to be certain that no saturated pixels were included in the final averaged 
kernel, and to ensure high S/N, the 200 mV SPR dataset was chosen for this analysis (as 
opposed to 0, 100 and 300 mV SPR data). 

Figure M below shows the 9x9 pixel kernel resulting from the data taken at 40V and 
200mV SPR. Despite the fact that the outer 9x9 annulus in the kernel below appears to 
be statistically indistinguishable for the errors, we note that all values in this annulus are 
systematically above zero indicating that the noise floor has not been reached. For this 
reason, the 9x9 annulus is included in the kernel. The values listed represent the percentage 
of the total signal for an average kernel, based on approximately 64,000 events. Measurement 
errors represent the standard deviation across all SPR events. 
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Fig. 2. IPC model for 3x3 pixel array. xtalk00F2.eps 
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Fig. 4. Charge diffusion example images, with Qpixci/Qtotai for each pixel. 
xtalk00F4.eps 
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Fig. 6. Image of pixels that have been reset using the SPR method for esti- 
mating IPC. xtalk00F6.eps 




Fig. 7. Illustration of background area used for correction. The shaded region 
of 2-pixel width is the annulus used for correction. xtalk00F7.eps 



15 





SFR: Fractionat Voltage in Center Pixel vs V^^, at 4W 
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4.B.2. Model 

As can be seen in the IPC data kernel in Figure [H IPC can spread significantly beyond 
the 3x3 grid depicted in Section I2.BI In fact, for each detector tested using SPR, the IPC 
distribution extended to at least a 9x9 grid. Additionally, each detector type (where type is 
determined by changing design and processing) was found to have a unique IPC distribution 
signature. Based on these findings, the IPC model is an averaged set of SPR data for each 
detector type, and given in terms of a (the average fractional signal in nearest neighboring 
pixels, per Section [2.61) . The use of a is in keeping with the standard modeling approach for 



IPC and provides a variable to capture the standard deviation of the measured and summed 
SPR kernels in the ^^Fe model (discussed in Section |4]E]). Variation of a is also needed to 



model the effects of varying levels of optical crosstalk on centroiding precision. 

In the example above (Figure [9]), a = (4.51 + 4.61 + 4.61 + 4.51)/4.0 = 4.56. Using this 
a, an IPC model (Figure fTOl) was developed to capture the signal-spreading signature of the 
BO generation. Each cell is a factor of a: 

where the central pixel, C = (1/a) J2 (of the grid of coefficients). 

4-0. IPC asymmetry findings 

Given the ROIC readout asymmetry in these detectors, we specifically set out to determine 
if a corresponding asymmetry exists in the IPC. The JMAPS test camera at JPL is designed 
to read out the full 4096x4096 pixel detector using 32 equally spaced outputs. This means 
that the detector is virtually divided into 32, 128x4096 pixel columns. The columns are read 
out alternately from the right and left, resulting in all odd columns being read out in one 
direction and even columns in the other direction. Figure [TT] shows this left /right, odd/even 
dichotomy. 

We conducted a column-by-column comparison of the IPC charge distribution for BO 
detectors. To perform this analysis the hits from all even and all odd columns were separately 
averaged. The first and last columns of the array were excluded from this analysis, since they 
suffer high dark current at the edges that could contaminate the results. A clear horizontal 
dichotomy resulting from the direction of the readout is visible in the odd vs even columns, 
as shown in Figure [121 

We performed this analysis for all of the SPR data taken for 0-300mV Vresct values and for 
all 4 substrate voltages. Table [2] summarizes our findings. In the readout direction (horizontal 
direction) the offset induced by the readout is approximately 5% depending on whether the 
column is odd or even. In addition, a vertical bias is evident in every column (both even 
and odd) of approximately 3.5%. The fifth and eighth columns of Table [2] show the residuals 
resulting when hits from both odd and even columns are averaged. We find the residuals to 
be low in the horizontal/read direction after averaging odd and even columns, but a 3.5% 



1^ 



residual is still present in the vertical direction. These findings will be used in the selection 
criteria of ^^Fe hits, as described in Section I4.DI 
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Fig. 10. IPC model based on SPR data. xtalkOOF10.ps 



20 



"Odd" Readout, 128 pixels wide 
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Fig. 12. SPR-generated IPC kernels from even (left) and odd (right) columns. 
xtalk00F13.eps 



22 



Table 2. Results of the odd/even column systematic offset test 



Vsubl (V) 


Vreset^ (V) 


x^ 


d 


^error^ 


jI 


J^ 


h 

y error^ 


10.00 


0.30 


1.08 


0.93 


-0.23% 


0.94 


0.94 


6.01% 


10.00 


0.20 


1.07 


0.94 


-0.29% 


0.97 


0.97 


3.27% 


10.00 


0.10 


1.06 


0.94 


-0.28% 


0.97 


0.97 


3.34% 


10.00 


0.00 


1.05 


0.96 


-0.39% 


0.96 


0.96 


3.52% 


20.00 


0.30 


1.08 


0.93 


-0.29% 


0.97 


0.97 


2.87% 


20.00 


0.20 


1.07 


0.94 


-0.28% 


0.97 


0.97 


3.16% 


20.00 


0.10 


1.06 


0.94 


-0.27% 


0.97 


0.97 


3.32% 


20.00 


0.00 


1.05 


0.96 


-0.42% 


0.97 


0.97 


3.39% 


30.00 


0.30 


1.08 


0.93 


-0.33% 


0.97 


0.97 


2.96% 


30.00 


0.20 


1.07 


0.94 


-0.26% 


0.97 


0.97 


3.27% 


30.00 


0.10 


1.06 


0.94 


-0.26% 


0.97 


0.97 


3.38% 


30.00 


0.00 


1.05 


0.95 


-0.34% 


0.97 


0.97 


3.43% 


40.00 


0.30 


1.08 


0.92 


-0.28% 


0.97 


0.97 


3.07% 


40.00 


0.20 


1.07 


0.93 


-0.21% 


0.97 


0.97 


3.32% 


40.00 


0.10 


1.06 


0.94 


-0.24% 


0.96 


0.97 


3.49% 


40.00 


0.00 


1.05 


0.95 


-0.30% 


0.96 


0.96 


3.53% 


Averages 








-0.29% 






3.46% 



"substrate voltage 

''SPR reset voltage 

'^Xright/xieft for all even columns 

"^Xright/xieft for all odd columns 

'^% residual error from hits in both columns 

■'^yupper/yiower for evCU ColumUS 
''yupper/yiower for odd ColumUS 

^% residual error from hits in both columns 



4.D. ^^Fe data acquisition and analysis 

^^Fe measurements were taken at the same voltage settings as the SPR measurements, using 
the JMAPS camera. For these measurements the source was exposed for up to 12 seconds. 
During the first two seconds, a mechanical shutter shielded the detector from the source. This 
was done to measure the offset from the zero point, due to reset, without contamination from 
the ^^Fe hits. In the remaining 10 seconds the shield was removed and the detector was non- 
destructively read out every 2 seconds, resulting in 5 frames of a progressively increasing 



23 



population of ^^Fe hits. 

The charge diffusion component of crosstalk is heavily dependent on the depth at which 
the Kq or K^ photon is absorbed. Two things must be taken into consideration when selecting 
good hits. First, the primary goal is to select hits that are as perpendicular to the detector 
surface plane as possible. Secondly, the hits must be selected so that all absorption depths 
are sampled. We define an accurate measurement of ^^Fe as the average over all absorption 
depths. To ensure this, we have used the following method to define normal ^^Fe hits: 

1. First the image was offset corrected. To avoid any noise associated with settling, the 
frames used for this analysis were the fourth and fifth frames. The fourth frame was 
taken as the offset frame and subtracted from the fifth frame. 

2. A small range of central pixel values were identified as hits. The range chosen is in 
ADU and strongly depends on the gain of the system. Therefore, for this initial rough 
estimate, the range was taken directly from a subsample of hits on the image without 
regard to the symmetry of the hit. The sum in a 9x9 region around the central pixel 
represents the majority of the analog digital units (ADUs) resulting from the hit. A 
Gaussian distribution centered on the Kq, peak results when these sums are plotted in 
a histogram. We defined the median of a Gaussian fit to this distribution divided by 
the total electrons expected for a Kq, hit (1620e~) as the gain of the system. In our 
case the gain was 1.19 +/- 0.04 ADU per electron. On average ~ 23, 000 hits proceed 
to Step 3 as potential K^ hits. 

3. Double hits or hits falling inside the parameter of the targeted hit region, can skew 
the crosstalk results. To remove double hits, all of the potential Kq, hits from Step 2 
were centered and stacked in the z direction, a process depicted in Figure [131 Kernels 
with double hits were removed by fitting a Gaussian distribution along the z-axis of 
the stack for each pixel except the central 3x3 in the kernel stack. Hits in the stack 
with pixel values greater than the mean +/- 3o" over the stack were considered double 
hits and the hit was removed from the analysis. The remaining hits in the stack are 
considered single hits and were used in the following steps of the analysis. 

4. A local dark correction was then performed on the remaining single hits. The pixels in 
an annulus of 2 pixels outside of the 9x9 central hit are defined as background pixels 
(see Figure [7]). The median of the pixels in this annulus was subtracted from the values 
of the inner 9x9 pixels. 

5. The symmetry criteria are perhaps the most challenging criteria to incorporate into the 
^^Fe reduction method. Noise and detector effects prevent even a perfectly symmetric 
hit from appearing symmetric in the data. This leads to confusion of off-center hits 
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with perfectly centered hits contaminated by noise. To choose perfectly centered hits 
we use the following criteria: 

(a) The noise of the system can be measured by fitting a Gaussian distribution to all 
pixels in the image except those affected by ^^Fe hits. The sigma of this distribu- 
tion represents the 1-sigma noise of the measurement. Two times this sigma is 
the maximum difference in counts that the left/right and upper /lower pixels can 
have, due to the noise. For all of our measurements, this noise was ~ 15 ADU 
(la). 

(b) As discussed in Section 4.3, the BO detectors show both horizontal and verti- 
cal systematic readout asymmetry. The asymmetry in the horizontal direction is 
column-dependent while the asymmetry in the vertical direction is not. We used 
the normalized kernels from our IPC analysis to estimate the offset in ADU ex- 
pected to result from the asymmetry induced by the readout. In most cases this 
asymmetry in effect contributed ~ 5 — 10 ADU to the total noise. 

(c) With the noise and asymmetry tabulated we defined centered hits as those hits 
having left/right and upper/lower differences that are less than 2-sigma plus the 
asymmetry offset. In most cases this resulted in a total permissible variation 
between left /right and upper /lower pixels of ~ 35 — 40 ADU. 

6. Finally, in order to ensure only K^ hits are included in the final crosstalk kernel (for 
the purposes of comparing the data to the model), the sum over all 9x9 pixels of each 
symmetric hit was taken. Using the gain from step 2, any hit with a sum from the 9x9 
that had values less than or equal to the K^ escape peak or greater than or equal to 
the K/s peak [10] were excluded from the stack. We note that the peaks of the K^ lines 
are used, rather than the edges of these peaks (based on Janesick [10]), in choosing 
our range. We do this because the noise is high enough for our system that the K/3 
peaks cannot be distinguished from the K^^ peak when plotted. To avoid preferentially 
excluding Kq, hits that overlap with the K^ hits due to noise, we established this 
definition. We recommend using only Kq, hits for those systems with noise that is low 
enough to distinguish the three peaks. 

7. After applying these criteria, the remaining hits were averaged to generate a 9x9 kernel 
(in the case of the 40V data) in units of ADU. The crosstalk kernel was defined as the 
charge in a single pixel divided by the total charge in the 9x9 kernel. 

8. This test was repeated 9 times, resulting in 9 kernels for all 9 datasets. The standard 
deviation across all nine measurements for each pixel represents the error on the final 
kernel due to statistical differences between each data set. 
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This method describes the reduction strategy for the 40V ^^Fe data. With lower substrate 
voltage settings we find that a kernel larger than 9x9 is needed to capture all of the escaped 
charge from the central pixel. We discuss this in more detail in Section 14. F[ 

4.E. ^^Fe 4OV data and model 

4.E.I. Data 

Figure [Ml shows the ^^Fe kernel resulting from the acquisition and analysis strategy described 
in Section 14. Dl for the H4RG-10 BO detector. Kernel pixel values are the percent measured 
signal for an average of 9 tests. The total number of good hits used for this kernel was 1094. 
The errors are also shown for each pixel value and represent the standard deviation of the 
pixel values over all 9 datasets. We define pixel crosstalk as the value of the average four 
nearest neighbor pixels (shown in grey) divided by the central value (shown in yellow). For 
the H4RG-10 BO detectors we obtain a measured crosstalk value of 9.83% for this ^^Fe data 
set. 

4.E.2. Model 

To generate the ^^Fe model, the average IPC kernel (/ipc,ij) , as described above in Section 
4.2, is used to represent the IPC spreading typical of the detector. The uncertainty associated 
with the standard deviation (o"cp) of the averaged single-pixel-reset (SPR) results is included 
in the ^^Fe model through convolution of the charge diffusion kernel with a set of IPC kernels 
derived from new alphas (anew)- These anew values are calculated from a set of central pixel 
values (Cncw) selected from a Gaussian distribution built around the averaged central pixel 
value (Cave) and the standard deviation of that averaged value, cXcp: 

Cnew = cTcp [-2 In (1 - rni)]^^^ cos {2v:rn2) + Cave (14) 

anew = ^^^ 7 (15) 

l^ij JlPC.ij 

where, rni and rn2 are random numbers in the interval (0, 1) (uniform deviates). The 
convolution is performed for a suitably large number of anew values, with noise added to 
each convolved ^^Fe kernel. (Noise values are randomly selected from a Gaussian distribution, 
with a noise standard deviation (o"noisc), 

anoise = Y^ (16) 

where G is the measured gain, e^ is the expected background noise (in e~), and 1620 is the 
number of charge pairs generated by a K^ photon absorbed in silicon.) The ^^Fe kernels are 
summed, and the final set of pixels is normalized so that the sum of charge across all pixels 
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in the final kernel = 1. The two-dimensional convolution of the charge diffusion and IPC 
kernels is performed using the Python routine convolve2d.py. Figure [15] shows the results of 
the ^^Fe model. 

A comparison of this ^^Fe model kernel with the measured ^^Fe data in Figure UM indicates 
that the combination our SPR model and charge diffusion model predicts the actual ^^Fe 
kernel very well at 40V. The central values for the model kernel are within the errors of 
the actual ^^Fe data as both tables above show. The SPR and charge diffusion model also 
predicts the values in the surrounding pixels very well. We note that some of the pixels 
around the central pixel do not quite match the model to within the errors. This is likely 
because the signal in the surrounding pixels is much less than the central pixel. When these 
pixels are normalized to the total in the kernel, the effect of noise in the normalization is 
much higher for these pixels than the central pixel. Thus, in addition to the statistical errors 
across all hits, there is error in the normalization due to the noise of the system. A system 
with lower noise than the camera presented here will likely find model values even closer to 
the observed values. We explore the predictive power of our model in more detail in the next 
section. 

4.F. ^^Fe modeling across voltages 

^^Fe measurements were taken at substrate voltages ranging from 10 to 40V. This section 
describes the crosstalk results for all voltage settings and compares these results with the 
model predictions. 

When extracting the kernel for the 10, 20, and 30 V measurements, it is critical to take into 
consideration broader spreading induced by charge diffusion. For 10, 20 and 30 V settings, 
we found that 15x15, 13x13, and 11x11 square pixel kernels respectively were needed to 
capture all of the charge from the ^^Fe hit. Because the gain of the system should not change 
with substrate voltage, the gain from the 40V measurement was used for the lower voltage 
settings, since here the Kq, distribution is much tighter. 

Table [3] shows the central and neighboring pixels resulting from the ^^Fe measurements 
and those predicted by the model. 
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Table 3. Fe model and data results for substrate voltages 10-40V 



Model and Data Comparison for 10-40V ^^Fe 


Central Pixel 


Vsub(V) 


kernel size 


model 


data 


std dev 


10 


15x15 


32.24% 


29.16% 


2.61% 


20 


13x13 


47.63% 


46.85% 


2.79% 


30 


11x11 


55.57% 


54.65% 


1.65% 


40 


9x9 


60.51% 


60.43% 


1.58% 






X 


-right Pixel 


10 


15x15 


9.52% 


8.88% 


0.57% 


20 


13x13 


7.79% 


7.65% 


0.38% 


30 


11x11 


6.48% 


6.63% 


0.27% 


40 


9x9 


5.72% 


5.92% 


0.17% 


X-left Pixel 


10 


15x15 


9.50% 


8.88% 


0.25% 


20 


13x13 


7.77% 


7.67% 


0.41% 


30 


11x11 


6.48% 


6.57% 


0.36% 


40 


9x9 


5.71% 


5.94% 


0.16% 
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Fig. 13. Double hits are removed by fitting a Gaussian to stacked pixels. Pixels 
with crosses were excluded from iteration. xtalk00F14.ps 
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0.00 
±0.06 


0.03 
±0.05 


0.03 
±0.06 


0.00 
±0.05 


0.06 
±0.06 


0.07 
±0.04 


0.00 
±0.06 


0.03 
±0.04 


0.02 
±0.05 


0.02 

±0.06 


0.01 
±0.06 


0.10 
±0.07 


0.07 
±0.07 


0.13 
±0.04 


0.07 
±0.08 


0.05 
±0.05 


0.06 
±0.06 


0.01 
±0.05 


-0.01 
±0.07 


0.03 
±0.07 


0.18 
±0.11 


0.39 
±0.05 


0.66 
±0.07 


0.35 
±0.07 


0.16 
±0.06 


0.05 
±0.06 


0.04 
±0.08 


0.01 
±0.07 


0.14 
±0.06 


0.36 
±0.09 


1.69 
±0.05 


5.76 
±0.20 


1.70 
±0.06 


0.39 
±0.06 


0.08 
±0.07 


0.03 
±0.08 


-0.03 
±0.11 


0.10 
±0.09 


0.66 
±0.06 


5.94 
±0.16 


60.43 
±1.58 


5.92 
±0.17 


0.69 
±0.10 


0.12 
±0.06 


0.04 
±0.10 


0.03 
±0.06 


0.16 
±0.06 


0.43 
±0.07 


1.89 
±0.05 


5.87 
±0.18 


1.79 
±0.04 


0.45 
±0.08 


0.10 
±0.09 


0.01 
±0.07 


0.03 
±0.06 


0.07 
±0.06 


0.16 
±0.08 


0.43 
±0.08 


0.69 
±0.10 


0.41 
±0.08 


0.16 
±0.08 


0.05 
±0.05 


-0.02 
±0.07 


0.03 
±0.06 


0.02 
±0.05 


0.06 
±0.06 


0.07 
±0.02 


0.12 
±0.07 


0.10 
±0.06 


0.03 
±0.05 


0.01 
±0.06 


-0.02 
±0.06 


0.01 
±0.07 


0.03 
±0.05 


0.02 
±0.05 


0.03 
±0.08 


0.05 
±0.05 


0.01 
±0.06 


0.03 
±0.05 


0.02 
±0.07 


0.02 
±0.06 
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0.01 


0.02 


0.04 


0.04 


0.04 


0.03 


0.03 


0.02 


0.01 


0.03 


0.04 


0.07 


0.12 


0.14 


0.13 


0.06 


0.04 


0.02 


0.02 


0.07 


0.17 


0.39 


0.66 


0.39 


0.19 


0.08 


0.04 


0.06 


0.09 


0.40 


1.67 


5.67 


1.64 


0.40 


0.12 


0.04 


0.03 


0.13 


0.61 


5.69 


60.51 


5.72 


0.60 


0.14 


0.03 


0.06 


0.15 


0.42 


1.73 


5.79 


1.75 


0.43 


0.13 


0.07 


0.05 


0.10 


0.18 


0.42 


0.69 


0.40 


0.19 


0.10 


0.02 


0.03 


0.04 


0.09 


0.10 


0.14 


0.12 


0.05 


0.03 


0.01 


0.00 


0.02 


0.03 


0.06 


0.05 


0.05 


0.04 


0.02 


0.02 
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Figure [T6] contains plots of the central pixel value and nearest neighbor value, both for the 
measurements and for the model. While the model comes close to the results we have for the 
lOV data, it does not predict the data to within the errors. At lOV the distributed charge 
inhabits a very large kernel: 15x15 pixels, for which the charge is very close to the noise floor 
of the system. We therefore conclude that the lOV data presented here is not high fldelity. 
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5. Summary 

While the design of source- follower CMOS-hybrid arrays eliminates some noise sources such 
as charge transfer inefficiency (CTI), read noise and dark current, it also results in an in- 
crease in IPC, in contrast to other detectors. IPC can be problematic for some astronomical 
applications like astrometry and photometry, which relies on high S/N to obtain accurate 
and scientifically useful results. 

In this paper we have used the SPR method to model IPC in TIS H4RG-10 arrays. We 
combined a charge diffusion model with this IPC model to understand the global effect of 
crosstalk on these arrays. Our primary findings are as follows: 

• The IPC in the H4RG-10 BO detectors extends out to 9x9 pixels based on our single 
pixel reset method, requiring a extension to Moores IPC model to capture all of the 
charge spread. 

• The single pixel reset method shows a read out bias that is column dependent. This 
effect introduces an inherent systematic bias in the IPC pattern. This bias must be 
included in the ^^Fe reduction strategy in order to extract all good hits from a ^^Fe 
image. 

• We have developed a model for combining IPC and charge diffusion effects as a function 
of substrate voltage. Our model results closely match our empirical results. 

• Our results suggest that the SPR method alone combined with the model of charge 
diffusion described here can provide enough information to simulate the dominant 
properties of H4RG-10 BO crosstalk to high accuracy without need for a radioactive 
source. In addition the SPR data with this charge diffusion model can be used to very 
accurately predict astrometric and photometric performance [4]. 

• We find that the H4RG-10 BO generation of detectors shows crosstalk of ~ 10%, 
and a signal loss of ~ 40% in the central pixel at 40V. A risk reduction effort is 
currently underway to reduce this crosstalk by half in a new generation (Bl) of H4RG- 
10 detectors. 

The model based on SPR data presented here reproduces the ^^Fe data to within the statis- 
tical errors and noise properties of our system. In a companion paper we use the models to 
explore the effects crosstalk at various levels (5-15%) have on astrometric and photometric 
applications [4]. 



34 



6. Acknowledgements 

The research described in this paper was performed in part by the Jet Propulsion Laboratory, 
California Institute of Technology, under contract with the National Aeronautics and Space 
Administration. The authors are very grateful to Chris Paine, Allan Runkle, Stuart Shaklan, 
and Larry Hovland, for all of their help designing, building and operating the JMAPS H4RG- 
10 astrometric camera at JPL. We are also very grateful to Andre Wong, Bill Tennant, and 
Yibin Bai at TIS for their informative thoughts and suggestions that have greatly improved 
this paper. 

References 

[I] J. Janesick, J. T. Andrews, and T. Elliot, Fundamental performance differences between 
CMOS and CCD imagers: Part 2, SPIE Proceedings, 6267 (2006). 

[2] B. N. Borland, et al.. Laboratory and sky testing results for the TIS H4RG-10 4k x 4k 
10-micron visible CMOS-hybrid detector, SPIE Proceedings Vol. 6690 (2007). 

[3] B. N. Borland, et al.. Initial laboratory and sky testing results for the second generation 
H4RG-10 4k X 4k, 10 micron visible CMOS-Hybrid detector, SPIE Proceedings Vol. 7439 
(2009). 

[4] J. Hubbs, et al.. High performance 4096x4096 pixel imaging sensors for the JMAPS 
Astrometry Mission, presented at the Military Science Symposium, Orlando FLorida 
USA 2011. 

[5] B. N.Borland, US Naval Observatory, 3450 Massachusetts Avenue, Washington, B.C. 
20392, USA, et al are preparing a manuscript to be called " Effects of pixel crosstalk on 
astronomical measurements using a large format CMOS-hybrid detector." 

[6] L. M. Simms, Hybrid CMOS SiPIN Detectors as Astronomical Imagers (ProQuest, UMI 
Bissertation Publishing, 2011). 

[7] L. Cheng, "Interpixel Capacitive Coupling," (thesis) RIT Bigital Media Library Repos- 
itory, https://ritdml.rit.edu/handle/1850/8414. 

[8] A. C. Moore, Z. Ninkov and W. J. Forrest, Interpixel capacitance in non-destructive 
focal plane arrays. Focal Plane Arrays for Space Telescopes, Proceedings of the SPIE, 
Vol. 5167, pp. 204-215 (2004). 

[9] S. M. Sze, Physics of Semiconductor Devices, Second Edition (Wiley-Interscience Publi- 
cation, 1981). 

[10] S. E. Holland, et al., Bevelopment of Back-Illuminated, FuUy-Bepleted CCB Image 
Sensors for Use in Astronomy and Astrophysics, in Proceedings of IEEE Workshop on 
Charge-Coupled Devices and Advanced Image Sensors , (1997), pp. R26-1 R26-4. 

[II] J. R. Janesick, Scientific Charge- Coupled Devices , (SPIE Press, 2001). 



35 



[12] http://henkle.lbl.gov/optical_constants/getdb2.html 



[12] S. Seshadri, D. M. Cole, B. R. Hancock and R. M. Smith, Mapping electrical crosstalk in 
pixelated sensor arrays, in High Energy, Optical, and Infrared Detectors for Astronomy 
III, Proceedings of the SPIE, Vol. 7021 , D. A. Dorn and A. D. Holland, ed. (2008), pp. 
1-10. 

[13] G. Finger, R. Dorn, M. Meyer, and L. Mehrgan,Interpixel capacitance in large format 
CMOS hybrid arrays, SPIE 6276, 62760F (2006) 



36 



IPC Kernel Mode! (% alpha) 




0.90 



' 0.75 



-0.60 



0.45 



0,30 



0.15 



"0.00 



0.1 - 

S£ 

^ 0.09 

■mO-0«5 - 

2: 

« 0.08 
^ 0.075 

e 

"1 t>.07 ■ 

'g 0.065 

C 

o 

t; o.(jft 

w 

u 

u. 

0.055 


Fraction of Signal in Nearest Neighbor for Fc^^ Model & Data 


















1 










\ 












) 


\\ 


V 














■ 


^ 


^■ 




















5^ 
















1 


N 


^ 


















N 


i 




















1 

i 














^**N..,„_""*"^i,„^^ T 




































0.05 1 1 1 1 1 ... ■ 
5 10 15 20 25 30 3S 40 45 

Subslrutti Voltagi; (V) 



